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Near  the  resonance  (uo  %  u>p)  a  high  amplitude  plasma  wave 

decays  into  a  purely  growing  electromagnetic  perturbation  and 

two  Langmuir  wave  sidebands.  The  instability  has  a  growth  rate 
2  2  2  2 

Y  -  ujp  kQlv0l  /8k  v*  for  a  Langmuir  pump  of  finite  wavenumber 

2  2 

kQ  >  k  vj c,  and  y  ^  u>p | VQ |  /16c  for  a  laser  light  pump.  In 

the  latter  case,  the  instability  has  an  inhomogeneity  threshold 
2  3 

| vq |  Wp  L/2  c  *  1  and  is  a  source  of  spontaneous  generation  of 
magnetic  field. 


; 
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I.  INTRODUCTION 


Oscillating  two  stream  instability,  in  which  a  plasma  wave  decays  into 
a  zero  frequency  electrostatic  mode  and  two  Langmuir  wave  sidebands  has  been 
known  for  some  time  and  shown  to  be  an  important  nonlinear  process  leading  to 
soliton  formation,  plasma  collapse,  etc.  In  this  paper  we  have  investigated 
the  electromagnetic  oscillating  two  stream  instability  of  electrostatic  and 
electromagnetic  plasma  waves  of  finite  wave  number  by  which  a  pump  wave  decays 
into  a  purely  growing  electromagnetic  mode  and  two  sidebands,  viz.,  Langmuir 
waves  of  shorter  wavelengths.  For  an  electrostatic  pump,  i.e.,  a  Langmuir 
wave  (ko[|^),  the  low  frequency  perturbation  has  a  finite  longitudinal 
component  of  electric  vector  (i.e.,  k  *  ?  #  0)  whereas  for  an  electromagnetic 
pump  (kQ  _L  or  a  dipole  pump,  the  perturbation  is  purely  electromagnetic 
(i.e,,  it  •  £  *  0).  In  the  latter  case,  it  is  similar  to  the  Weibel  insta¬ 
bility  in  an  anisotropic  plasma  and  leads  to  spontaneous  generation  of  magnetic 
field.  Bodner  and  Eddleman^  had  some  time  ago  pointed  out  this  instability 
of  intense  laser  radiation  by  considering  a  dipole  pump  and  a  homogeneous 
plasma.  Here  we  have  studied  the  electromagnetic  oscillating  two  stream 
instability  of  Langmuir  waves  of  finite  wavenumber;  kQ  t  0  introduces  major 
qualitative  and  quantitative  effects  on  the  parametric  instability.  In  the 
case  of  an  electromagnetic  pump,  i.e.,  a  laser  radiation,  we  have  studied  this 
instability  including  the  effect  of  plasma  inhomogeneity  and  shown  its  im¬ 
portance  as  a  source  of  self-generated  magnetic  fields  in  laser-fusion 
experiments.  ^ 

Self-generated  magnetic  fields  of  a  few  megagauss  have  been  observed  in 

2—6 

a  number  of  experiments  and  computer  simulations  of  laser  produced  plasmas. 

The  origin  of  these  fields  is  believed  to  be  due  to:  (1)  resonance  absorption, 
in  which  the  density  oscillations  due  to  the  p-polarized  light  couple  with 


the  drift  velocity  to  produce  a  time  average  current  and  hence,  the  magnetic 

9-12 

field;  (2)  ponderomotive  force,  which  contains  a  curl  free  term  due 

to  collisions;  (3)  mutually  perpendicular  gradients  in  density  and  temperature 

13 

7n  x  VT;  (4)  thermal  instability  and  (5)  Weibel  instability  of  counter 
streaming  electron  beams.  For  high  power  density  lasers  the  electromagnetic 
oscillating  two  stream  instability  should  become  an  important  source  of 
magnetic  field  generation. 

The  mechanism  of  the  instability  is  as  follows.  A  high  amplitude  plasma 
wave  (o>o#  &o)  propagating  through  a  plasma  produces  an  oscillating  drift 
velocity  of  electrons  vq  in  the  direction  of  the  electric  field  We 

perturb  the  system  by  an  electromagnetic  perturbation  (w,  lc);  (lc  ]_  £q)  and 
the  magnetic  vector  5  is  perpendicular  to  both  lc  and  2  .  For  an  electro¬ 
magnetic  pump  £  |  |lc,  where  as  for  an  electrostatic  pump.  In  the 

former  case  the  magnetic  field  of  the  perturbation  couples  with  the  oscillatory 
drift  velocity  of  electrons  to  produce  a  vq  x  ?  force  driving  the  high 
frequency  sidebands,  viz.,  electrostatic  plasma  waves.  The  sidebands  couple 
with  the  pump  to  produce  a  low  frequency  current  driving  the  instability. 

For  an  electrostatic  pump,  the  low  frequency  perturbation  has  finite  density 
fluctuations  also,  which  lead  to  an  important  nonlinearity  producing  the  side¬ 
bands.  The  sidebands  interact  with  the  pump  to  produce  a  low  frequency 
ponderomotive  force  as  well  as  a  transverse  current  driving  the  instability. 

In  both  the  cases  the  sidebands  are  coupled  through  the  relativistic  velocity 
dependence  of  mass  and  hence  the  coupling  coefficient  for  the  parametric 
decay  is  modified. 

In  Section  II  we  have  studied  the  electromagnetic  oscillating  two  stream 
instability  of  Langmuir  waves  of  finite  wavenumber.  In  Section  III,  the 
instability  analysis  is  extended  to  an  electromagnetic  pump,  viz.,  laser 
radiation,  including  the  effect  of  a  density  gradient.  A  discussion  of  results 


is  given  in  Section  IV, 
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II.  INSTABILITY  OF  LANGMUIR  WAVES 


We  consider  the  propagation  of  a  Langmuir  wave  x  Eq  exp-i(o)ot  -  kQ  x) 

in  an  uniform  plasma.  This  produces  an  oscillatory  velocity  vq  =  e  io/m  i  o>0> 

a  density  oscillation  n  =  n°  e  k  E  /m  i  w2  t  and  decays  into  a  low  frequency 
J  ooooo 

electromagnetic  perturbation  t  =  c  t  x  l£/u>  and  two  high  frequency 

electrostatic  sidebands  $ .  (w .  . ) ;  j  =  1,2,  co1  9  =  w  -  w  ,  It  9  =  k  -  k  , 

it | | z  and  k  >>  kQ.  The  pump  and  sidebands  produce  high  frequency  current 

densities,  J_2x^w  ~  ^o^  =  ~ni  e  v0^  an<*  +  ^Wo^  "  ““n2  e  Vo^  also> 

2  2  2 

which  produce  a  magnetic  field  B  2y  =  ^  ^  J_2x/C(k  ~  ),  important 

for  deriving  the  sidebands. 

The  response  at  (w,k)  and  2*^1  2^  *s  °bta*nec*  by  solving  the  equations 
of  continuity  and  motion. 


H  +  V  •  (nv)  =  0,  (1) 

■—(mv)  +  v  •  V  (mv)  =  -el£  -  -vxS--Vn  ,  (2) 

3 1  c  n  - 

m  =  m  (1  +  v  •  v/2c2) , 
o 

We  expand  the  quantities  v  and  n  with  subscripts  o,  1,  2  referring  to  the 
pump  and  two  sidebands  and  the  quantities  without  a  subscript  referring  to 
(u>,1c)  component.  For  u>^,  k^  component  Eq.  (1)  gives 

n-  -  -k  v  n/2u +  n°  k  v.  /o>  .  (3) 

loo  1  o  lz  1 

v,  is  obtained  from  Eq.  (2)  retaining  the  relativistic  effects  through 
lz 

the  first  term  on  the  left-hand-side  and  v  x  I  nonlinearity  through  the 
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second  term  on  the  right.  The  last  term  in  Eq.  (2)  can  be  expressed  in  terms 

2 

of  by  using  the  Poisson’s  equation  k  <t>^  =  -4Tren^.  Solving  these  coupled 
equations  along  with  the  similar  equations  for  the  upper  side  band  we 

obtain 


v  0  k  v  0 

.  '  o'  ,  ,  2ir  T1oo  27rne 

“ - v  J  + - r~ 

II  0  2  ^  ,  12  o  x  m  4  2 


8c 


c  ok 


(4) 


where 


I*  r 


e2^2  '  '  „  2  V1  "  ,  2 2  "x  "  a) 


2nv  k  v  _ 

o  _  o  o  2irne 

vj  “ 


8C 


k*-  * 

c 


o  k 


£  =  1  —  — ^  — 
j  2  2 

J  U).  (l) 

J  O 


e  +  ° 


(5) 


,  2  ,22  |  ,2  ,2,,  2,2 

co  kv  |  v  1  k  +  4  u>  /c 

o 


Q  2  .2  ,  2,  2  9 

8c  k  *■“  4  oo  /  c 

o 


j  =  1,2,  to2  =  4TTn°e2/in  (1  +  3|v  |2/8c2), 


In  writing  these  equations  we  have  expressed  B  *  y  4tt1  J^/kc.  In  Eqs.  (4) 
and  (5)  the  first  terms  on  RHS  are  due  to  relativistic  mass  correction,  second 
terms  due  to  vq  x  B  force  and  the  last  terms  due  to  nvQ  nonlinearity  in  the 
equation  of  continuity.  The  last  term  vanishes  for  a  dipole  pump  (k^  «  0) . 


For  U),k  component  Eqs.  (1)  and  (2)  give 


eE 

x 

/  -  — —  , 
x  n  k 
o 


(6) 


kv  eE  .  . 

- +  — t  +  ir~  (-k  v  v1  +  k  v  v0  ) , 

a)  o  m  ioo  2co  o  o  lz  o  o  2z  ’ 
n  o 
o 


n  -  - 


eE  n°  k  v  e ($..  -  <t>~)n0 
z  o  o  o  l  2  o 


m  ikv^ 
o  e 


2v 


mu 


(7) 


(8) 


Here  we  have  Ignored  the  relativistic  corrections  as  <o,k  is  not  a  normal 


mode.  Using  the  Poisson's  equation  and  Eq.  (8),  one  obtains 


6 


o  , 

n  e  k  v 

o  o  o  , , 

n  *  ~T~  ~2^~  ^1 

O 


2 

0)  J 

^  ~2~ 
o>  e 


where 


U) 


2 
U) 

^Ei 
.2  2  2 
k  v  a) 
e 


(9) 


(10) 


The  low  frequency  current  density  may  be  written  as 


T  o 

J  =  -n  e  v 

X  ox 


nn  e  v  /2  -  n.  e  v  /2  * 
1  o  2  o 


o  2 
n  e  E 
o  x 

mio> 


k2v 

TT<*1 


(ID 


Using  Eq.  (11)  in  the  wave  equation. 


V2E  +  SL.  E  *  -  Auiw  j 
x  2  x  2  X 

c  c 


one  obtains, 
E 


_ia>_  V*!  ~  ^2} 

“  2c2  1  +  u)2/k2c2 

P 

,!AV. 

x  "  Q  .  2 2  2N  * 

8tt(1  +  w  /k  c  ) 

P 


(12) 


(13) 


Employing  Eq,  (13)  in  Eqs.  (4)  and  (5),  we  obtain  the  nonlinear  dispersion 
relation 


ef  e* 
le2 


[v  ty/l6c^ 
o 


(14) 


*  -  r 


.  2  2  2  „ 
-  .  k  c  u)  2 

_ I  +  _2 _ pii 

.  .  2. .2  2  2  ,  2  2  2  J 

1  +  h>  /k  c  k  v  w  e 

p  e 


+  A 

i  2 


ei,2  “  1  - 

“P 

2 

“1,2 

.2  2 
k  v. 

K 

a  -  e 

j.  1  < 

°  2  . 

[- 


,2  ^  ,  2/  2  ,  2  2  2  0 

k  +  4w  /c  k  c  w  ,  2 

_ 9 - +  _2 _ Eli 

2,  2/,  2  2  2  J 


2  ,  2 1.  2  2  0  *,2  *  *  /  ^ 
1+u)  /k  c  2(k  -4uj/c)  kv  we 

p  o  e 


,1/2 


Defining  6  ■  w  **  u>o(l  +  4)  ,  Eq.  (14)  takes  the  form 
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2 

0) 


(15) 


The  largest  growth  occurs  for  6 
1 2 


0  i.e.,  for  A  ~  0  because  id  J  id  or 

o  p 


,2  2 

k  v 

e 


U) 


k 


k2  +  4cj2/c2 
o 


4c 


i  ,  2;i  2  2 

1  +  a)  /k  c 
P 


,22  2 

k  c  a) 

+  - Ei] 

2  ,  2,  2.  .2  2  2  J 

2(k  -  4w  /c  )  k  v  0)  e 

o  e 


2u> 

(A) 

O 


(16) 


2  2  2  2 

Eq.  (15)  can  be  solved  in  different  limits.  For  k  c  /k  v  <1,  the  last 

o  e 

terms  in  the  expressions  of  i|>  and  A  can  be  neglected  and  the  growth  rate 
y  =  -loo  turns  out  to  be 


Y  =  “p^ol  2/16c;2 


(17) 


the  same  expression  as  obtained  by  Bodner  and  Eddleman  in  the  limit  of  a 

dipole  pump.  The  low  frequency  mode  is  purely  electromagnetic  in  the  limit 

of  k  -*  0. 
o 

2  2  2  2 

In  the  opposite  limit  of  kQc  /k  >  1,  i.e.,  a  pump  wave  of  substantial 
wavenumber,  the  last  terms  in  ij;  and  A  are  predominant;  these  terms  appear 
through  the  density  fluctuations  associated  with  the  low  frequency  mode,  viz., 
the  last  terms  in  Eqs.  (4)  and  (5).  On  further  assuming  a>  <  Zw^m/m^  (m^  is 
the  mass  of  ions)  the  dispersion  relation  yields, 


Y 


o 

8 


I  2  2 

v  k 
1  o  o 

2,2 
v  k 
e 


(18) 


It  must  be  mentioned  here  that  the  low  frequency  mode  is  a  mixed  mode  with 

2  1/2 

magnetic  field  B  -  2  w  .  E  /t o(ww  )  >  E  >  E  .  For  u>  >  2  oj  m/m. ,  the 

°  y  pi  z  p7  z  x  p  i 

dispersion  relation  becomes  a  bicubic  equation. 
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III.  INSTABILITY  OF  LASER  RADIATION: 

A  SOURCE  OF  SELF-GENERATED  MAGNETIC  FIELDS 


Here  we  consider  the  colinear  propagation  of  laser  light  pump 
x  Eq  exp-i^^t  -  k^z)  and  all  the  decay  waves  along  the  z  axis,  i.e., 
ko|ik|jz„  The  z  components  of  the  low  frequency  ponderomo t ive  force  and 
the  current  density  and  hence,  the  density  oscillations  vanish,  i.e.,  the  low 
frequency  mode  is  purely  electromagnetic.  The  nonlinearity  driving  the  low 

frequency  perturbation  comes  from  J  =  -n^ev^/2  -  n^ev^/2  [the  same  as  Eq.  (11)] 

->■ 

and  the  sidebands  are  driven  by  the  vq  x  B  force.  Sideband  potentials  are 

given  by  Eqs.  (4)  and  (5)  with  n  =  0.  The  nonlinear  dispersion  relation, 

for  k  >>  k^,  turns  out  to  be  the  same  as  Eq.  (14)  with  4;  and  A  having  kQ  =  0. 

2  2 

The  uniform  medium  growth  rate  is  then  given  by  Eq.  (17),  viz.,  y  =  |  v^  |  /16c  . 

11  “1  15  2  2  1/2  2 
For  a  Ndtglass  laser  y  -  3  x  10  sec  at  4  x  10  (1  -  w  /o>  )  W/cm  ,  i.e., 

P  o 

the  growth  time  is  ^3p  sec.  It  may  be  worthwhile  mentioning  here  that  the 

growth  rate  of  electrostatic  oscillating  two  stream  instability  is  given  by 

222422  2  22  22 

Y  =  Imu),  a)  -  k  c  -  gj  / (u)  -  oj  . )  =  2w  .  S  6/ (gj  -  6  )  where  3  “  k  |v  |  /4w  . 

1  s  pi  pi  o  o  1  o'  p 

At  large  power  densities,  the  growth  rate  maximizes  for  6  =  y  and  y  = 

max  max 

2  2  2  2 
k  Anco  |v  |  /8v  . 

D  p  1  o  '  e 

Now  we  consider  the  effect  of  a  density  gradient  on  this  instability. 

All  the  waves  are  taken  to  propagate  along  the  density  gradient  i.e.,  z  axis. 

kc  9  9 

Using  Eqs.  (4)  and  (5)  and  writing  Jx  =  •£— r  By,  E1  =  -  —  4^,  E2  =  -  —  <t>2 
one  obtains 


-S  Ei +  <- 

dz 

d2E„  J  - 

~T  +  (_L“ 

dz  v 


2  2  ..22 
U>,  “  0)  v 

_ E  _  -  SL ...  o 

2  2  2  hl 

v  v  8c 


e 

2 


e 

,2 


,2  2  2 

J  1  u)  v  (ij 

£ 1  _£_  _  o  o  p 

2  2  z  ~  2c  2  y 

v  8c  v 

e  e 

.2  2 
w 


72>  r2 

8c 


o 1  o 

v2  8c 
e 


V  o> 

2  Ei  =  2?  “2  By 
v  J 
e 


(20) 


(21) 


fate*-1 
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Eq .  (11),  in  compliance  with  the  Faraday's  and  Ampere's  laws,  gives 

,2  2  0 

d  B  u)  v  2 

_ Z  -  P  R  = _ 9.  (v  -  v  \ 

2  2  By  2c  ,  2  (EI  V* 

dz  c  dz 


(22) 


Eqs.  (20)  -  (22)  are  very  similar  to  those  obtained  for  oscillating  two  stream 

instability  suggesting  that  the  electromagnetic  instability  is  absolute.1^  To 

calculate  the  threshold  for  the  instability,  we  take  u>  =  0  (the  growth  vanishes 

jv  | 2 

at  the  threshold)  and  assume  co  =  (to  — u)^)(l  +  Z/L).  Then  Eqs.  (20)  -  (22) 

p  o  ,  _  z  o 


4c 


take  the  form 


d2S  ,  Z 
dz" 

D 


v 

1  c 
4c  \ 

D 


)  S  =  0 


(23) 


d"  n 

Z  ^ 

0 

2  D  ~ 

.2  D  ~ 

9 

dz 

lXd 

2 

cAd 

,  2 

O) 

V 

b 

- §■  B  ® 

O 

dz2  y 

c2  y 

2c 

(24) 


(25) 


dz 


2  2  2 

where  X_  ■  v  /w  ,  S  *  E,  +  E„,  D  =  E,  -  E„.  Eqs.  (23)  -  (25)  can  be  solved 
D  e  p  i  l  l  i 

employing  the  Fourier  transform  technique  (cf.  Liu^) 


et  -  7- 
1  4  7T 


dk[exp-i(kz  F  Lk^X^/3)] 


f  lVol 

<  a?  exp(i  L  - k)  +  a, 

l  4c  1 


r.  _  lVol  ^  -1  kc.  \ 

exp[i  L  - (k - tan  — )]  r 

o  2  c  w  I 

2c  p  ; 


(26) 


E2  =  }  dk-IexP“i(Icz  +  Lk^A^/3)] 

”00 

/  lv0|2 

*  I  a^  exp(i  L - k)  -  a^ 


l 
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V  u  1  \, 

f.  ,  o'  p  t  -1  kc, 

expli  h - 7r  (k  — -  tan  — J 

„  ^  c  w 

2<*  p 


(27) 


At  large  values  of  -z(-z  >  L*v  j^/2c^)  the  integrals  can  be  evaluated  by  the 
saddle  point  method  with  saddle  points  at  k  -  +  (-z/L)^^A  \  corresponding 
to  ingoing  and  outgoing  waves.  Since  the  source  of  these  waves  is  situated 
at  z  -  0,  there  cannot  be  any  Incoming  wave  at  z  *>  This  demands^ ^ 


v 

o 


2  0)  L/2c  ! 
P 


(2n  +  1)  ;  n  is  an  integer. 


14  2 

The  corresponding  threshold  power  density  of  a  Nd: glass  laser  is  ^5  x  10  W/cm 

for  a  density  scale  length  of  100  pm.  Here  we  have  taken  into  account  the 

enhancement  of  the  pump  field  (which  goes  as  an  Airy  function)  in  the  vicinity 

1/4 

of  resonance,  E  ^ (L/A  )  *  E0  (vacuum). 

As  the  amplitude  of  the  instability  grows  with  time,  the  magnetic  field  B 
attains  large  values  and  modifies  the  dispersion  relation  for  plasma  waves, 
increasing  the  fn  euency  shift  6, 

6  =  (w2  +  w2)1/2  -  a)  (1  +  A) 1/2 
pc  o  ' 


where  wc  =  eB/mc.  The  instability  would  be  completely  stabilised  when 
2  2  _ 

6  a,  I  v0  I  /16c  >  i.e.,  'v  ojp  |  VQ  j /2/2c .  This  corresponds  tc  Gauss 

15  2 

at  laser  intensity  "vlO  W/cm  .  However,  there  may  be  other  saturating 
mechanisms  also  for  tit  is  instability. 
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IV.  DISCUSSION 


Electromagnetic  oscillating  two  stream  instability  constitutes  an 

efficient  channel  of  decay  of  Langmuir  waves  to  shorter  wavelengths  and  may 

play  an  important  role  in  the  saturation  of  beam-plasma  instability.  The 

decay  waves  propagate  almost  perpendicular  to  the  pump  wave  and  hence,  one 

should  consider  the  effects  of  finite  beam  size  on  this  instability.  The 

effect  of  finite  kQ,  wave  number  of  the  pump,  is  to  excite  density  fluctuations 

along  with  purely  growing  electromagnetic  perturbation.  These  fluctuations 

provide  strong  interaction  between  the  sidebands  and  the  pump  and  thereby 

2  2  2  2 

enhance  the  growth  rate  y;  y  %  w  k  I v  J  /8v  k  for  k  >  k  v  /c. 

p  o  1  o 1  e  o  e 

For  a  high  power  laser  radiation,  this  instability  is  an  important  source 

of  spontaneous  generation  of  magnetic  fields,  which  have  been  observed  in  a 

number  of  experiments.  The  effect  of  finite  kQ  is  negligible  and  the  growth 
2  2 

rate  y  -  w  |vq|  /8c  .  For  a  Nd: glass  laser,  the  growth  turns  out  to  be  a  few 
15  2 

pico  second  at  MO  W/cm  .  The  density  threshold  for  this  instability  is 
14  2 

M  x  10  W/cm  for  a  density  scale  length  MOO  tm.  As  the  amplitude  of  the 

instability  (i.e.,  of  the  magnetic  field)  grows,  the  four  wave  resonance 

is  detuned  and  the  instability  is  stabilized  at  a  level  cj  ^  w  v  /2/2c.  This 

c  p  o 

corresponds  to  the  magnetic  fields  of  a  few  megagauss  at  laser  intensities 

15  2 

MO  W/cm  . 
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